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SUMMARY 

The paper presents a mathematical discussion of the 
laminar boundary layer, which was developed with a view 
of facilitating the investigation of those boundary layers 
in particular for which the phenomenon of separation 
occurs. The treatment starts with a slight modification of 
the form of the boundary layer equation first published by 
von Mises. Two approximate solutions of this equation 
are found, one of which is exact at the outer edge of the 
boundary layer while the other is exact at the wall. The 
final solution is obtained by joining these two solutions 
at the inflection points of the velocity profiles. The final 
solution is given in terms of a series of universal functions 
for a fairly broad class of potential velocity distributions 
outside of the boundary layer. Detailed calculations of 
the boundary layer characteristics are worked out for the 
case in which the potential velocity is a linear function of 
the distance from the upstream stagnation point. Finally 
the complete separation point characteristics are deter- 
mined for the boundary layer associated with a potential 
velocity distribution made up of two linear functions of 
the distance from the stagnation point. It appears that 
extensions of the detailed calculations to more complex 
potential flows can be fairly easily carried out by using 
the explicit formulae given in the paper. 

1. INTRODUCTION 

The theory of the laminar boundary layer has two 
important applications: First, the computation of skin 
friction in the range of low Reynolds Numbers; second, 
the explanation of the separation phenomena. Both 
problems require the calculation of the "development 
of the boundary layer" under the assumption of a 
given pressure distribution along the wall. In some 
simple cases, as in the case of constant pressure and in 
that of a pressure decreasing proportionally to a given 
power of the length measured along the wall, the partial 
differential equation of the boundary layer theory can 
be reduced to a total equation and integrated without 
difficulty. In these cases aU "cross sections" of the 
boundary layer show "similar velocity profiles" and 
essentiaUy only the scale of the velocity and the length- 
scale over the cross sections, i.e., the "thickness of the 
boundary layer" are variable. However, in general, 



especially in the case in which the pressure is increasing 
along the wall, the distortion in the shape of the veloc- 
ity profile is the very point of interest. In these cases 
the solution of the partial differential equation itself is 
necessary. Unfortunately, all methods indicated until 
the present involve such an amount of numerical work 
as to discourage any engineer anxious to use the theory 
for solution of practical problems. 

K. Pohlhausen (reference 1), at the suggestion of the 
first author, tried to reduce the boundary layer prob.- 
lem to the solution of a total differential equation in the 
following way: He chooses a certain family of plausible 
"velocity profiles" with variable shape in such way 
that the influence of the rate of change of the pressure 
along the wall is taken into account by a boundary 
condition, but one parameter (for instance, the "thick- 
ness" of the boundary layer) is left undetermined. 
Then he applies the momentum law to a strip of the 
boundary layer enclosed between two adjacent cross 
sections using the integral relation introduced by the 
first author. This relation leads to a total differential 
equation for the undetermined parameter and by 
integration of this equation the "development of the 
boundary layer" is obtained. 

This method has been criticized especially by von 
Mises (reference 2), who pointed out that the results 
depend greatly on the choice of the family of velocity 
profiles. This is true to some extent, in spite of the 
fact that in the simple cases mentioned above plausibly 
chosen velocity profiles lead to olose agreement with 
the exact solutior both for the value of the skin 
friction and the thickness of the boundary layer. 
Furthermore, the fact that the method can be applied 
to the case of the turbulent layer, in which we ignore 
the partial differential equations of the motion, but 
know the approximate shape of the velocity profile, 
has to be considered as a great advantage of the 
method. As a matter of fact, all more or less successful 
developments in the theory of the turbulent boundary 
layer are based on the Karman-Pohlhausen procedure. 

However, in the case of the laminar boundary layer, 
the equations of the motion are known and the diffi- 
culties are merely mathematical. Thus a method which 
contains less arbitrary assumptions than the Pohl- 
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hausen procedure and involves less numerical calcula- 
tions, than the methods used by Hiemenz, Boltze, 
Thorn, Green (reference 3), etc., would represent a 
desirable improvement of the situation. 

The authors became interested in the solution of the 
boundary-layer equation in connection with their in- 
vestigations concerning the maximum lift (stalling 
point) of airfoils. It is generally known that stalling 
is due to the separation of the boundary layer at the 
upper surface of an airfoil, but in general it is assumed 
that the boundary layer is turbulent, except in the 
immediate neighborhood of the stagnation point, over 
the range of high Reynolds Numbers corresponding to 
actual flight or usual wind-tunnel test conditions. 
Recent investigations in the 10-foot diameter wind 
tunnel of the Guggenheim Aeronautics Laboratory at 
the California Institute of Technology have revealed a 
systematic and very considerable influence of the 
amount of turbulence in the wind stream on the value 
of the maximum lif t (reference 4) . This influence could 
be correlated very exactly with the influence of the 




Figure I.— Fluid element in the boundary layer. 

turbulence on the so-called "critical" Reynolds Number 
of a sphere (Reynolds Number characteristic for the 
sudden drop of the drag coefficient). As L. Prandtl 
(reference 5) pointed out, the influence of turbulence 
on the latter phenomenon is due to the fact that the 
regime in the boundary layer changes earlier from 
the laminar to the turbulent state, if the turbulence of 
the wind stream is increased. The authors concluded 
that a similar process might explain the influence of 
turbulence on the stalling of airfoils. The conception 
may be stated as follows: The low values of lift maxima 
in quiet air streams are due to the separation of the 
laminar boundary layer, while in a turbulent stream 
the regime changes before separation occurs and so the 
lift maximum is raised. For the investigation of such 
an hypothesis, the calculation of the development of 
the laminar boundary layer from the stagnation point 
up to the separation point appeared as a necessity. 

This paper presents the method suggested by the 
authors for the computation of the development of a 
boundary layer, especially suitable for the case of in- 
creasing pressure involving separation. As will be 
seen, the fundamental equation of motion is used in the 



form given by von Mises in Ms paper (reference 2) 
with a slight modification. This modification furnishes 
the possibility that the solution of an equation identical 
with that of ordinary heat conduction, furnishes a 
first approximation for the velocity field in the region 
away from the wall. In addition, the shape of the 
velocity profile near the wall is calculated from the 
given pressure distribution along the wall. The joining 
of these two solutions is carried out in a rather rough 
way, at least from a mathematical point of view. The 
authors hope that experts in applied mathematics will 
accomplish the task of improving their somewhat 
rudimentary procedure. 

2. THE FUNDAMENTAL EQUATION IN THE THEORY OF 
THE BOUNDARY LAYER 

All of the following investigations are restricted to 
the case of 2-dimensional stationary motion. We 
use a system of curvilinear coordinates. The lines 
^= constant represent stream lines, so that the value 
of the parameter f gives the amount of fluid flowing 
through between the wall (^=0) and the particular 
streamline corresponding to the value of $ considered. 
A system of curves perpendicular to the stream lines 
■will be called the 7i-curves; the length measured along 
the ^-curves may be denoted by s, the curvature of 
the stream lines by kf, that of the 7i-curves by Je„. 

We introduce the following notation: 

•wis the magnitude of the velocity, the 
direction of the velocity being given 
by the direction of the stream lines. 

p is the pressure. 
<r s , <r„, r are viscous stresses acting on soctions 
normal and parallel to the stream 
lines (as shown in fig. 1). 

p is density of the fluid. 

v is kinematic viscosity of the fluid ==ju/p. 

Let us consider the equation of momentum applied 
to the component of the momentum in the direction 
of the stream lines: 

In this equation b/bn denotes the differential quotient 
along the 7i-curves, where the letter n is used for the 
length measured along the w-curve's. 

Let us consider the components of the viscous stress 
<r s , <T n , t. The normal stresses <r s and o-„ are propor- 
tional to the rate of extension of the volume element in 
the s and n directions, the shearing stress is propor- 
tional to the rate of shear. We write 



<r n =-2A« 
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Now the theory of the boundary layer is based on 
two fundamental assumptions: 

(a) that inside of the boundary layer the differ- 
ential quotients perpendicular to the direction of 

flow increase proportionally to -n=> as the kinematic 

viscosity v tends to zero, while the differential quo- 
tients in the direction of the flow remain finite. 

(b) that the curvature of the stream lines k$ remains 
finite if j»-»0. 

According to these assumptions we can neglect the 
terms containing <r s and a n in comparison with the 
terms containing t, we can neglect k^r in comparison 

with ^> and finally we can replace t by/x ^ 

Hence we obtain the relatively simple equation 



bs\'2 + p) v bn 2 



(1) 



In this equation p appears as variable both along 
the stream lines and the n-curves. Following Prandtl's 
method, we should investigate the variation of p 
in the direction perpendicular to the flow by con- 
sidering the equation of momentum applied to the 
component taken in that direction. According to 
the assumptions (a) and (6), we easily obtain the 
result that the variation of p perpendicular to the 
streamlines is of the order of -yfv. Hence in equation 
(1) we replace p by its value p 0 corresponding to 
71=00, i. e., outside of the boundary layer. In this 
way p appears as a given function of s. 

We may now consider somewhat more exactly the 
meaning of the two coordinates s and n used in equa- 
tion (1). s was defined as the length measured along 
the stream lines. However, it is obvious that with 
the approximation used in this theory we can use s as 
the length measured along the wall, and n the length 
measured along the normal curves (perpendicular to 
the stream lines) or along a straight line perpendicular 
to the wall. The approximation used in the theory 
does not permit of a discrimination between such 
different definitions of the "cross section of the 
boundary layer" (reference 2). 

We therefore replace the coordinates s and n by two 
parameters <p and f chosen in the following way: 

(a) Every boundary layer calculation starts from a 
given "outside flow" in which the influence of viscosity 
can be neglected, the influence of viscosity being 
restricted to the boundary layer itself. The velocity 
of the outside flow at the wall shall be given by U(s). 
We introduce the parameter <p by the relation 

<P = J Q U(s)ds, i.e. <p is the line integral of the "outside 

velocity" along the wall. If the outside motion is a 
potential motion, <p is simply the value of the potential 

along the wall. Obviously we have to write ^= U-^- 



(b) We replace n by the parameter ^ as defined 



above and write ^— = 
bn 



b$ 



Let us calculate 



bn 2 ' 



b 2 u b ( bu\ b ( b%\ b 2 fv?\ 

Hence we write equation (1) in the following form: 
TT b/u 2 \, TT ldp 0 b 2 /u 2 \ 

Now in the "outside flow" the Bernoulli equation 

d /p U 2 \ 
holds, so that p +~2~ ) = 0 and we obtain: 

T7 b{v?-U 2 \ b 2 fu 2 \ 

or talcing into account the fact that U is a function 
only of <p and independent of ^ 

b/ TP-% 2 \ u b 2 f U 2 -v? \ 
&A 2 ; v Ub^\ 2 / 

TP— ii? 

We call the quantity — ^ — =z the "energy defect." 

It can be interpreted as the loss of energy of a fluid 
particle of unit mass flowing along a particular stream 
line as compared with the energy of a particle outside 
of the boundary layer. Introducing z in the last 
equation, we obtain 



bz u b 2 z 
b v ~ v Ubif 



(2) 



We notice that this equation is nearly identical with 
the equation given by von Mises (reference 2). The 
difference lies essentially in the use of the parameter 
<p instead of the length s. It will be shown below that 
this modification facilitates the practical use of the 
equation. 

First it is obvious that the use of <p instead of s does 
not introduce any complication in the set-up of the 
problem, because in most cases equally simple expres- 



sions are available for the functional relations 



TP 



and for U=g(s). Second, it is easy to give to equation 
(2) the following interpretation: In the system <p, 4>, 
TP—u 2 

the energy defect z= — ^ — obeys an equation analo- 
gous to the equation of heat conduction. The only 
coefficient appearing in the equation is the quantity 

lb 

v jj> i.e., the kinematic viscosity reduced in the ratio 

jj> which is the ratio between the local velocity inside 

the boundary layer and the outside velocity at the same 
cross section. 
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Neglecting first the influence of the "reduction fac- 
tor " jj> i.e., replacing v jj by v, we obtain a first approx- 
imation for the solution of the fundamental equation 
(2). Because the approximation is most accurate in 

11/ 

the outside portion of the boundary layer, where -jj 

is nearly equal to unity, we will refer to this solution 
as the "outside solution" and denote it with z a . It 
has the great advantage that it can be obtained by 
using the well-known methods for integration of the 
heat-conduction equation. 

3. DIMENSIONLESS FORM FOR THE FUNDAMENTAL 
EQUATION 

Before proceeding to the calculation of this first 
approximation, a dimensionless form will be given for 
"the fundamental equation (2). 

We introduce the following notation: 

TJ 0 is a suitably chosen characteristic velocity 

involved in the problem. 
L is a suitably chosen characteristic length 
involved in the problem. 

R = is the corresponding Reynolds Number. 

Then <?*, u*, TJ*, z*, shall denote dimensionless 
quantities defined by the following notation: 



* v_ 

* 'HZ 



* U 0 L2 



u 

u u 0 



(3) 



Substituting these values in equation (2), we obtain 



dz* 



1 u* d 2 2* 



"4 U* b^* 2 

where can be expressed by z* in the following form: 

u* _ I 2z* 
U* y 1 u* 2 

For convenience, we omit the stars in the succeeding 
analysis, so that the simple letters denote dimension- 
less quantities. The same letters are therefore used 
in the future in a different sense from that in which 
they were used in the previous sections. 

The final equations can now be written: 



dz_ lu b 2 z 
d*> _ 4 Ubi/, 2 



(4) 



and 



u__ 

U~\ U 2 



(5) 



It is necessary to indicate the boundary conditions 
for the equation (4). Obviously the value ^ = 0 corre- 
sponds to the wall and ^= eo to the transition to the 

U 2 

outside flow. Therefore, «=0, or z=-^- = z 0 (say) 

for ^=0 and 2=0 for <» . It is found convenient 

U 2 

in the future to use z 0 and interchangeably. 

4. GENERAL EQUATIONS FOR THE FIRST 
APPROXIMATION 



(OUTSIDE SOLUTION) 

Let us consider the equation 

02, 



■6) _ 1 C) 2 2ft) 

bcp 4 OiA 2 



(6) 



which we obtain from (4) replacing by 1, or neglect- 

2z 

ing jp in (5) in comparison with 1, and writing z — z a . 

This approximation cannot in general hold near the 

2z 

wall, where jp approaches the value 1. However, it 

is not without interest to discuss this simple approxi- 
mate solution, which we shall later use for the outside 
part of the boundary layer. 

Considering z as "temperature", and replacing <p by 
the time t, $ by a reduced length-coordinate x, the 
following analogous heat-conduction problem can be 
formulated: The temperature at the end point x = 0 
of a bar of infinite length is a given function of the 
time t beginning at t = 0. At infinite distance (x = °° ) 
the temperature is equal to zero. The temperature 
distribution as a function of time and distance from 
x=0 is to be computed. 

For either heat conduction or boundary layer state- 
ment of the problem, we have to find the solution of 
(6), which satisfies the boundary conditions 2 U = 0 

for f= oo and 2 m = 2q= ^ ^ f° r 1^=0. This solution 
is given by the definite integral 

which is given in any text on the theory of heat conduc- 
tion. Hence with this approximation the velocity 
distribution in the boundary layer can be computed by 
quadratures, as TJ 2 is a given function of the outside 
potential <p. 

It is interesting that a relatively simple explicit 
condition can be found for the location of the separa- 
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tion point. Separation occurs if the slope of the 
velocity at the wall changes sign. Thus the condition 



bz 



for separation can be written in the form ^; = 0 for 
^ = 0. We obtain by differentiation 



For ^=0 the second term vanishes, and the first 
term becomes indeterminate, but its value can be 
determined by a partial integration. In this way we 
obtain 



fbz u \ = 2_ 



I 



dfc 



Thus the separation criterion is given as 
where 

r dC7 



» ^ dg 



d£=0. 



(8) 



Assuming ^ = 0 as stagnation point (Z7=0) and 
assuming positive values for £7, we easily see that 

df7 

separation can only occur after changes sign, i.e., 

after the velocity U has reached a maximum value 
and is decreasing. 

We will find later that this separation condition is 
not of great practical value, because as we approach 
the separation point, the approximate solution becomes 
very inexact near the wall. 

5. THE VELOCITY PROFILE AND DIFFERENTIAL 
EQUATION IN TERMS OF w=w (x, y) 

The differential equation (6) determines the energy 
U 2 —u 2 

defect z= — ^ — as function of the parameters <p 

and ^. For many applications, expressions for the 
velocity profile are needed, i.e., one has to determine 
the velocity u as function of the coordinates x and y, 
where x is the distance of an arbitrary cross section 
from the starting point of the boundary layer, and y 
is the distance from the wall. Instead of the real 
lengths x and y, dimensionless quantities x* and y* 
may be defined in terms of the characteristic length L, 
putting 

.y 

'L 



x*= T y* : 



For convenience we use the letters x and y instead of 
the letters x* and y*, as we did above for the other 
quantities involved in this investigation. Then it is 
obvious from the definitions given in equation (3) 
that u, x, y can be expressed by the formulas: 



u=-y/U 2 -2z=U- } Jl~ 
Jo U-Jo 

2 r*difr 2 p 



z_ 

Zo ) 



(9) 



These formulas enable us to calculate the velocity 
profiles corresponding to any given solution z(<p, 

It is interesting to transform the approximate equa- 
tion (6) into an equation for u, containing x and y as 
independent variables. Obviously the symbol bjb<p 
means in this case the rate of change of any quantity 
with increasing x while keeping constant y, while b/b<p 
means differentiation by keeping ^ invariable. There- 
fore 



_b_ 

dp 



D\&> das *yf 



where ^ is the inclination of the line ^= constant in 

the x, y coordinate system. The inclination ^ can 
be expressed by 

fog v 
btp~ il 
by 

where v is the (dimensionless) velocity in the y direc- 
tion. We obtain in this way 

bz a _ 1 f bz w , v bz a \ _dU 1/ bu bu\ 
b<p~U\bx + ubyJ~ dx U\ w bx +V by/ 

On the other hand, 

b 2 z b_( bu\ -2 b fbu\ -4 b 2 u 

by ~ bA u a VS a A*u) ?>y 2 

(from (9)). 

The equation (6) is therefore transformed into 



bu , bu 1 Ub 2 u , TT dU 
bx by Jiu by 1 



Ax 



Calling the dimensional quantities Ui, v u U u x u 
and yi for the moment, we obtain 



It! -5- — = v ^ — 9+ U\-j — 

0x1 by 1 u x byi axi 



(10) 



Comparing equation (10) with the t( exact" boundary 
layer equation deduced first by L. Prandtl: 

bui , bUi d 2 tti, rT dE7i 

bxi by! byf dxj v 1 

we notice that equation (10) exaggerates the influence 
of the term corresponding to the viscous friction, 
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multiplying the term v -^-f with the ratio — • We 
especially notice that for y=0, i.e., at the wall, u=0 



and so 



U b 2 u 



can only have finite value if 



u by 2 ulujr "~ vo ¥ ^ UD " by 2 



0- 



Hence the solution of equation (6), called in this paper 
the "outer solution", is composed of velocity profiles 

with vanishing second differential quotients at the 

wall. This is the reason that it cannot be expected 
that S shape velocity profiles, such as occur in a flow 
against increasing pressure before separation, will be 
satisfactorily represented by the solution of the 
approximate equation (6). 

6. POWER SERIES EXPANSION OF U* (*>) 

The analytical evaluation of the integral in (7) for 
an arbitrary function U 2 (<p) is in general difficult, if 
not impossible. Numerical or graphical methods are 
likewise of little use. We therefore adopt the pro- 
cedure of approximating to the function U 2 (<p) by a 
polynomial of the form 



i=0 



(12) 



With such a function the indicated integrations can be 
analytically carried out. In proceeding with the cal- 
culations to this end, it is convenient first to transform 



the expression for z u in (7) by introducing a new 
variable of integration defined by 



vV-£ 

A little, calculation gives at once 



Zu) — 



-JttJ±_ 

-fp 



(13) 



or introducing the particular function U 2 given in (12): 

z «-is s 'X v "(*-F)' d <'- 

Expanding the term in parenthesis by the binomial 
theorem and integrating term by term, we obtain: 



where 



2 C°e~ p ' 
Jrfr) = i= -55rd/3, and 0.' = 1. 



By repeated partial integrations J T where r is any in- 
teger greater than zero, may be expressed in terms of 
J 0 . Carrying out the partial integrations and col- 
lecting terms, we obtain, finally: 



Jq(x) 



z, 

where 
0o (a) 



9l (*) = 7= xe ^ + 3? Jo ( x ) 



J 0 (x) = 1 — %- \ e-? d/3, and 0 ! = 1. 



(14) 



The universal functions g r (x) can very readily be 
calculated using the standard tables of the probability 
integral appearing in J 0 (x). Once determined they 
are independent of the velocity function TJ(tp), pro- 
vided only that the latter has the polynomial form 
(12). Equation (14) therefore, in connection with 
(9), gives the complete analytical solution for our 
first approximation to the boundary layer problem 
for any case in which the potential velocity is such 
that £7* may be expressed as a polynomial in <p. 

Substituting (12) in the separation criterion (8), 
integrating term by term, and again applying succes- 
sive partial integrations, it is very easy to show that 
the separation criterion takes the form, where <p s is the 
value of <p at the separation point: 



J=l 



(2i-l)l' 



0 



(15) 



Hence the determination of the separation point is 
reduced to the solution of an algebraic equation with 
given numerical coefficients. 

7. DOUBLE POWER SERIES EXPANSION OF U 2 (*>) 

In principle we can approximate with any desired 
exactness to any practically important velocity func- 
tion TJ 2 (ip) by a polynomial of the form (12). For 
many important cases, however, such as, for example, 
the potential velocity around an airfoil, the number of 
terms required in the polynomial for a satisfactory 
approximation becomes very large. Since the ex- 
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pression for z u becomes rapidly more and more com- 
plicated as the number of terms in the polynomial 
increases, it follows that the solution of our problem 
becomes very laborious and complex in such cases. 
For this reason a solution has been investigated for the 
case in which U 2 (<p) is approximated by two distinct 
polynomials of the form (12) instead of by one. This 
solution, which is discussed in the present section, is 
somewhat more complicated formally than is (14), but 
is vastly simpler for carrying out actual calculations. 

We approximate to U 2 (v) by two polynomials, one 
of which we use for <p^<pi and one for where <p x 

is any fixed value of <j>, i.e.: 



U 2 (<p) =M<P) =SM for 



i=0 



U 2 (<p) =M<p) for v^ei 

i=Q 



(16) 



The procedure is indicated schematically in figure 2, in 
which the solid line represents the U 2 (<p) actually used. 

"We substitute (16) in the original integral for z a , 
change the variable of integration, expand, integrate 
term by term, and perform successive partial integra- 
tions just as in the analysis leading to (14). The 
final result may be written in the form: 



where 



(17) 



Zt = z 0} of formula (14) 

Si,* = &offo( ; ^ ) 
\-\<p— <Pi/ 



zp* is identical with z b * except b t — and n — >ra. 
0!=1 and g r (x) is as given in equation (14). If <p^<pi 
we use the simpler formula (14). 
Proceeding as before, the separation criterion now 
takes the form 



w . . ^"-'ilfr-l)! , 

(2i-l)l »"» + 



+ 



/ Vl ' 2 2 "i(2i 



2r)\ 
r)W 2 



0!=1, the second sum is to be taken to the 
larger of n or m, and b t =0 for \>n, >3 t =0 for £>m. , 



(18) 



For many practically important cases it is sufficiently 
accurate to approximate to U 2 (<p) by two power series 
involving no terms of degree higher than 3 in <p, i.e., 

hi=b s = =/34 = ft= =0. Insuch 

cases the formula (17) can be written in the following 
form, which is very convenient for purposes of 
calculation: 



z M (<p, $ = {b 0 (h 0 -h 0 *) + Mo*+7i<Pi£ri* 

-Y2¥>i. 2 fir2* + Y3¥>i 3 03*} 
+ (h - A,*) + PA* + 2y 2 <Pi (*7x* 

+ 02*)-373?l 2 (<72*+<73*)} 
+ 2fif 2 * + ff3*)} 

+ <P s {h(h 3 -h 3 *) + P3h 3 *} (for vzzvj [ (19) 

where 

+3^2+^3, h*=go*, h*=ga*+g*, etc. 
If tp^pi put Pi=gi*=hi* = 0 

Before proceeding with the general case, it is inter- 
esting to compare this approximation with the exact 




Figure 2.— Double power series representation of TJ l (s>). 

solution given by Blasius for uniform external flow 
along a flat plate. 

8. APPLICATION TO BLASIUS' CASE U= constant 

For this simple case, we have in accordance with (16) 

U 2 = 60, bi = b 2 — =j3i=0, <p± — co 

Hence from (19) and (14) 

zM) = boh = b 0 g 0 = |°[l - ^= ^e-P d/}] 

If we take U= Z7 0 =the characteristic velocity used 
in defining the Reynolds Number, then U has the 
numerical value 1, i.e.: 

For convenience, we denote the probability integral 
by P, i.e.: 

P( X ) = A C* e -p d p 
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Then from the basic equations (9), we have for the 
dimensionless velocity and distance from the wall 



vV 



2v; rv? 



dj8 



■yjR Jo *JP(fi 



(20) 



We arbitrarily define the boundary layer thickness 5 
as that value of y at which the velocity has 99.5 percent 
of the value which it has in the outside potential flow, 
i.e., 

y=8~u= 0.995 (21) 
From the tables of the probability integral P, we find 
that this occurs for -t=-= 1.8. 



Hence 



_2jh r- 

jRJo 
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Figube 3.— Dimensionless Telocity profiles for flow along a flat plate. 

The evaluation of the integral in (20) has been carried 
out analytically for 0^ -%=■ ^ 0.1 using the standard 



power series expansion: 



PW 



3 + 10 + 



•)• 



For 



0.1 the integration has been accomplished 



graphically with a Coradi integraph. Before present- 
ing the results it will be convenient to introduce the 
distance along the plate in a direction parallel to U 
and measured from the leading edge, wbich is also 
taken as the origin of <p. In order to fix ideas, we 
consider a plate whose total length in the direction 
of U is taken as the characteristic length L which 
occurs in the Eeynolds Number. Then in order to 
make our distance along the plate a dimensionless 
variable, we divide it by L. This dimensionless dis- 
tance we denote by a; in accordance with customary 
practice, so that O^a;^ 1. For the present case 



<p= Ux=x since 27=1. Hence <p may be replaced 
throughout by x. 

The results for u=u(y) are given in figure 3 by the 
curve labeled" "first approximation". The Blasius 
profile is also plotted for comparison. The expression 

5 53 

for 5 from this first approximation is s = ~^ V^> 

5.5 

while the usual Blasius value is S Biasius — -yfx. 

Here the dimensionless length 5 is, of course, also 
defined in terms of the standard length L. The shear- 
ing stress at the wall t 0 , as obtained from (^) wnll nas 
the following expressions: 

/ r 0 \ 0.282 / T 0 \ _ 0.332 

1st approx. Blasius ■yjRx 

It has been pointed out above that our first approxi- 
mation corresponds to an exaggeration of the viscous 
forces near the wall. The method of linearizing the 
differential equation introduced by Oseen retains the 
exact viscous terms but approximates to the inertia 
terms. It is of some interest, therefore, to compare 
the results given by the two methods. This is especi- 
ally easy for the present case of flow without pressure 
drop, since the Oseen method leads to an equation of 
exactly the same nature as does our first approxima- 
tion. The Oseen solution has, therefore, been worked 
out and will be briefly discussed here. 1 

- We start Avith the dimensionless form of the Prandtl 
equation which, for 27=constant, becomes (cf. equa- 
tion (11)): 

bu , bu b 2 u 
We then introduce the velocity defect ^ defined by: 

1t2= U— u 

Substituting this in the Prandtl equation, the Oseen 
method of linearizing is based on the neglect of all 
second-degree terms in or v. Carrying out this 
procedure, we obtain 

bit? _ 1 b 2 U2 

bx~UR by 2 ' 

1.6.; tl heat-conduction equation of exactly the type as 
was deduced for our first approximation (6). The 
boundary conditions are: vk = U= 1 at y= 0 and ih — 0 
at y= oo. The solution corresponding to these bound- 
ary conditions is obtainable at once in terms of the 
probability integral as 

fBU y 



2 

or in terms of the original variable u (since Z7= 1) 

' After the present paper was completed, the authors discovered that tho Oseon 
method results presented in the remainder of this section had boon investigated 
from a point of view rather different from theirs and presented in a paper by N. A. 
V. Piercy and H. F. Winny "Tho Skin Friction of Flat Plates to Osoen's Approx 
imation", Proc. Boy. Soc. Lond., Series A, Volume 140 (1933), p. 613. 
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. p/V-g y \ 

This velocity profile has also been plotted in figure 3 
and labeled "Oseen." The boundary-layer thickness 
defined as above is given by 

while the velocity gradient at the wall can readily be 
shown to be just twice that given by our first approxi- 
mation, so that 

/ to \ = 0.564 

\pU 2 /osem -JRx 

Comparing all of these results, we see that our first 
approximation gives excellent agreement with the 
Blasius exact solution as far as 5 is concerned, and a 
fairly satisfactory agreement with respect to velocity 
profile and wall shearing stress. The Oseen method 
of linearizing the differential equation leads, in the 
present case, to much greater errors in all of these 
characteristics. 

9. FURTHER APPROXIMATION TO THE OUTER 
SOLUTION 

In order to improve the accuracy of our solution, it 
is natural to look for some method of successive approx- 
imation or iteration starting with the first approximate 
solution discussed above. Von Mises (reference 2) 
has suggested a method based essentially on successive 
iterations in the <p direction. If we use subscripts 1 
to denote our first approximation and subscripts 2 for 
a second approximation, then in view of the form of 
the exact equations (4), (5), von Mises suggests 
essentially that we write 

This procedure has been investigated for the Blasius 
case and found to lead to inconsistencies with the 
given boundary conditions of the problem. The diffi- 
culties seemed to be basic ones, so that this method of 
approximation was abandoned and the following alter- 
native procedure was considered: 
The exact differential equation may be written in 

the form 

bs m_1^2/ /i_22_-A 
b v 4drp~4c)^ 2 VV U 2 V 

Then the first approximation is obtained by neglect- 
ing the right side of the equation, i.e.: 

5s M 1 b 2 g M n 
4 b^ _U - 

The following second approximation immediately 
suggests itself: 

bz 2 1 b 2 2 2 1 b 2 z a f L 2z a \ , , n 

bi~l W^l W\S 1 ~TP- 1 ) =i ^> *>> sa ^- 



Then ^ is the solution of a generalized, nonhomo- 
geneous, heat-conduction equation, in which the right 
side fi, is a given function of <p, From any standard 
work on mathematical analysis 1 ), the formal solution 
of this equation satisfying the given boundary con- 
ditions can readily be found to be 

if" c v r ^-">' 

z 2 (<p, £) = z a (ip, — 7= drj \ d£L e 
-yjirjo Jo 

Unfortunately, the authors have found the evalua- 
tion of this integral form of solution too intractable 
for the method to be useful for the case of any general 
function U{ip). Only in the Blasius case, where a 
separation of variables is possible so that the partial 
differential equation reduces to an ordinary differen- 
tial equation with independent variable -7= , has a 

solution of this second approximation been calculated. 
Since the method is not pursued further in the present 
paper, the analysis leading to this solution is not 
presented. The velocity profile so deduced is, how- 
ever, included in figure 3, with the label "second 
approximation". The approximation to the exact 
solution is seen to be quite good. 

In spite of this initial success the entire method of 
successive approximations was abandoned for the rea- 
sons mentioned above. It appears, however, that it 
might furnish an interesting and perhaps profitable 
field of study for the mathematician. 

10. THE INNER SOLUTION 

It has already been pointed out that the approximate 
form of the boundary-layer equation which leads to the 
outer solution discussed above is exact at the outer 
edge of the boundary layer, but is in general far from 
correct at the wall. A natural procedure is therefore 
to look for another approximation to the exact equa- 
tion which shall be correct- at the wall but possibly 
inaccurate far from the wall. The solution of this 
equation might naturally be called the "inner solu- 
tion", and the final method of approximating to the 
exact and complete problem would then be to join the 
outer and inner solutions together in some way. 

A first approximation to this inner solution is very 
easily obtained, but it is convenient to first change the 
notation slightly in setting up the equation. We intro- 
duce a new independent variable defined by 

u 2 TP 
f =20-2=2" where 2 0 =-2"- ^ 

Then the basic differential equation (4) may be written 
in the form 

1 cf. for example, Goursat's " Corns d' Analyse HI, " eh. XXIX. 
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where 



20 



Using a subscript ( )< to denote the inner solution, we 
approximate to (23) by writing 

W Vr*V z a '\b<p) )\ ■ (24) 
A(,<p) = -4z 0 '-Jz~ 0 ] 

where J is a function of <p, ip, i; t chosen in such a 

way that (24) approaches (23) as the distance from the 
wall goes to zero, i.e., as i^-*0. 

Since at the wall we have f =|^=0, it appears that 

to get a first approximation for our inner solution we 

may take = 0- This solution, denoted by fr, will 

then be exact at the wall. Discussing this simple 
approximation, we have 

We consider this as an ordinary differential equation 
and integrate for <p= constant. It is convenient to 
introduce temporarily the auxiliary dependent variable 

W= (^jtff) > an< i to consider £ i, as independent variable. 
Then we obtain at once 



(25) 



or integrating, 
Now writing 



ldW = A 



and returning to the original variables f i x and 
A httle calculation gives finally, assuming 

r, I (^=o)=o. 

^^Ib^-(B^-2A^)^BJ+4A^]. (26) 

If Biz^O, we may invert this expression obtaining the 
following expansion about ^ = 0: 

f Ji -£ 1 *+!^|riP/»+ (5^0). (27) 

At the separation point (j^) ~ 0 an< i hence Bf 
— 0. We have exactly, therefore: 



ur or 



h=(i) 4/3 ^H 



separation point : (28) 
J5 t = 0 



The different analytical form for at and away 
from the separation point for ^->0 is very interesting. 
In view of the exactness of the equation for £<„ at 
^ = 0, any further approximations for must exhibit 
the same behavior near ^=0 as is shown in (27) and 
(28). 

11. THE IMPROVED INNER SOLUTION AND ITS 
JOINING TO THE OUTER SOLUTION 

In developing further approximations to the inner 
solution, we must be guided by the procedure which 
will be used in joining the two solutions together. A 
discussion of this procedure requires first a brief re- 
considering of the approximate outer solution. The 
exact and approximate outer equations are, respec- 
tively: 

d Vi _£.Y /2 ™a bz * ld2g " 

The two equations are identical at the outer edge of 
the boundary layer where z—0 and would at first 
glance appear to become more and more different as 
the wall is approached, since z in general increases 
continuously from 0 to z 0 . Fortunately, however, the 

two equations are again identical when ■^ = 0, i.e., 

at an inflection point of any s(^) profile. For the 
Blasius case of no pressure drop this inflection point 
occurs just at the wall. For external flows with pres- 
sure drop in the direction of U it is almost certain 
that there is no inflection point at all. For flows with 
pressure increase in going downstream the inflection 
point moves out from the wall so that the s(^) profiles 
develop an S-shape. In view of the nature of the 
two equations, it therefore appears that the outer 
solution furnishes, in most cases, a satisfactoiy ap- 
proximation to the exact solution for the region be- 
tween the outer edge of the boundary layer and the 
inflection points of the zty) profiles. Proceeding from 
these inflection points to the wall the accuracy of the 
outer solution apparently becomes rapidly more and 
more unsatisfactory. The procedure adopted in view 
of this situation is the following: For flows with pres- 
sure decrease or constant pressure the outer solution 
is considered as a satisfactoiy approximation to tho 
exact solution. For flows with pressure increase the 
outer solution is used from the outside of the boundary 
layer to the curve connecting the inflection points of 
the outer solution z u (ip) profiles. For the region be- 
tween this curve and the wall the partial differential 
equation is replaced by a family of ordinary differ- 
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enfcial equations in f. The inner solution is composed 
of the functions £ t (f) obtained by the solution of 
these equations and joined to the appropriate 
profile of the outer solution for each value of <p. In 
the following discussion of the inner solution and its 
joining to the outer solution, we are therefore con- 
sidering only flows with pressure increase in the direc- 
tion of the external flow. 
In joining the two solutions it seems necessary to 

require that -^p be continuous, in view of the major 

importance of this second derivative in the equations. 
We must, therefore, join the two solutions in such a 
way that the inflection points of the outer solution 
coincide with inflection points of the inner solution. 
Referring to equation (25), we see that for a given 



value of <p, 



d* 8 



has always the same sign for all 



values of In other words, the first approximate 
inner solution has no inflection points and is, therefore, 
quite useless for our present purpose. Some 15 differ- 
ent procedures for finding a satisfactory inner solu- 
tion and joining it to the outer one were investigated 
in all, the more important of which will be briefly 
discussed below. 

It is convenient to call the point at which the inner 
and outer solutions are joined the "joining point" and 
to use a subscript j to denote quantities taken at this 
point. It is also convenient to introduce an alterna- 
tive independent variable ! = -t=> and to consider £ 

vV 

and z sometimes as functions of <p and £ rather than 
of <p and ^. The variable if appears continually 
throughout the en the analysis and seems to be a natural 
physical variable of the boundary-layer problem. Its 

close relationship to the variable -H=- which gives 

->Jx 

the Blasius similarity law is quite evident. It should 
be noted that inflection points in the f, \f/ plane 
correspond identically to inflection points in the f , £ 
plane, so that our preceding discussion is valid in 
either plane. 

The first procedure investigated was an iteration one, 
i.e., f was taken from the first approximation 



as 



dip 



This had the effect of moving the inflection 



point in from infinity, but in the cases considered still 
gave fi,->z 0 , i.e., the inflection point still lay outside 
of the boundary layer. This method was, therefore, 
abandoned and in all the succeeding investigations 
the inflection or joining point was forced to occur in 
the neighborhood of the wall through the choice of 

(^~)' At the wall we must have (^£) =0 since 
f(^=0) = 0 for all £>'s, while at the joining point we 



1 /df\* 

see from equation (24) that -pijhp) =1 - Hence 

a series of simple assumptions for which sat- 

isfied these two conditions was discussed. The various 
attempts will be more clear if we restate the immediate 
problem of the inner solution a little more explicitly. 
We have a second order differential equation (24) for 
f i as function of ^ (or £), <p being considered constant. 
We should like to impose the following boundary con- 
ditions on the solution of this equation: 

At the wall, ^ = 0, ^=0. 

At the joining point, ^ = ^ (determined 
from the outer solution). 



w(ID,-(i); 



r. 

In addition, we must choose so that at the joining 

point (^Sji) = 0- We have here three boundary condi- 
tions for a second-order differential equation and must, 
therefore, abandon one. The condition at the wall is 
of essential importance, especially for the separation 
phenomenon, and must be retained. Of the three re- 
maining possibilities two have been discussed. First, 
the condition ? { j=t aj was abandoned. Second, the 
joining point was considered to be dete rmin ed by the 
condition that ti } =!; Uj and the value of was not 
required to be the same as that for the outer solution, 
i.e., in general, ipt^fay The discarding of the condi- 
tion on was n °t considered, largely because the 
second method was felt to give satisfactory results. 
When the first method was used (^~) was assumed to 

to be given by (^) = |~" -^ or tne second method 
two assumptions with regard to (^p) were consid- 

the more natural, but it unfortunately leads to the 
appearance of elliptic integrals which make the suc- 
ceeding calculations very awkward. The second is very 
fortunate in that the inner solution is readily obtained 
in explicit form and involving only elementary func- 
tions. 

In order to discuss these various alternatives a fam- 
ily of simple external flows of the form U 2 = b 0 + bi<? 
was used. Inaccuracies in the different procedures 
would be expected to reveal themselves most strikingly 
at the separation point. Accordingly, the separation 
point location and the corresponding velocity profile 
were calculated for all three of the cases mentioned in 
the preceding paragraph, and all three procedures were 
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found to give very similar results. The calculation, of 
separation point locations followed methods analogous 
to that described in the next section. In view of this 
lack of sensitivity of the separation point character- 
istics to variations in the method of joining, the second 
assumption of the second method was adopted, i.e. 

Before giving the analysis and results for this pro- 
cedure we shall briefly outline the essential elements of 
the method finally adopted. For a given value of <p 

we determine f„ A (jjjjfy and &y from tne outer solu_ 




r-f;-Origin for £, 
M" " " ii 
Figure 4. — Method of joining inner and outer solutions to obtain the complete 
joined solution. 

tion, using the condition (^f) = 0- ^ e tnen assume 
■^r(sr^) =-»/— and deterniine the inner solution 
using the boundary conditions 

r,-o at *=o, Qi) r (^) j •* r«-r-»-*> 

We then determine f fj corresponding to and define a 
quantity e by the relation 

For |"<£V we use the inner solution as a function of 
<P and £*• For f >f j we use the outer solution, but for 
a given value of f 0 we use, not ? w or f a but £ M +e or 
lAo+eV^- Our final approximate solution then has 
the form indicated schematically in figure 4 by the 
solid line. This procedure implies that the outer 
solution is one which does not correspond to the 
actual wall but to a fictitious one displaced out into 
the fluid by an amount determined by e. This dis- 



placement is a function of the distance downstream 
from the stagnation point, and in general, has its 
largest value at the separation point. The error 
introduced through the use of this displaced outer 
solution will naturally depend on the relative size of e, 
increasing as e increases. In the 13 velocity profiles 
TP(<p) which have been investigated using this proce- 
dure the value of £5 (corresponding to the boundary- 
layer thickness) at the separation point always lay 
between 1.5 and 2.0 while the largest value of e was 
0.026. Hence the displacement was at most of the 
order of 1.5 percent of the width of the boundary- 
layer profile. It is, therefore, believed that the error 
introduced by the displacement of the origin of the 
outer solution is of negligible importance. In this 
connection one further procedure which was investi- 
gated should be mentioned. In this a slightly altered 
value of z 0 was taken for the outer solution and the 
alteration was taken to be such that e vanished. This 
implies that a solution was obtained for a shghtly 
distorted external potential flow, but the solution was 
logically consistent. The results with this method 
were practically indistinguishable from those obtained 
with the e procedure, so that this more complex 
method was not carried further. 

12. ANALYSIS AND GENERAL RESULTS FOR THE FINAL 
JOINED SOLUTION 

We first derive the expression for the inner solution 

finally used. Substituting (^) = ^/j^ ™ e< l ua - 
tion (24), we have 

Integrating in the same manner as before, we obtain: 



where, as before, 

Integrating -once again and introducing the bound- 
ary condition f i (i/'=0) = 0, we have 



sm 



2A-JT} 



2.4 VfJ 



It is convenient to rewrite these expressions in terms 
of the variable £. Replacing A by its original expres- 
sion from (24) a little calculation gives for the final 
inner solution expressions: 
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+ sin" 



go - 4a 



where 



£ 2 =- 



(29) 



fj \<>yJv=o 



(30) 



We obtain the following expressions, giving the inner 
solution joining-point characteristics, by putting f!=^ 
in (29): 



+ sin" 



Vi+# 




(3D 



To join to the outer solution for any given value of <p, 
we first determine 



,— »and\-<nr/ from the outer 
solution (17) or (19) using the relation y-^" 0 / 

ALA / d £A 

= £ ; \"w )r \~w J i The second e ^ uation ° f (31) 

then determines B and this, with the preceding equa- 
tion, gives Finally, we obtain e from 

. (32) 

We now have the complete and joined outer and inner 
solutions for a given <p. The analysis may then be 
repeated for as many values of <p as are desired. The 
fact that the inner solution is given in the form £= £(f ,), 
and that inversion to f 4 (£) is very difficult, is of no 
practical importance, since what we finally want is the 
velocity profile u=u(y) for an arbitrary value of x, i.e., 
of ip. For considering only the inner solution for the 
moment, we may write the basic relations (9) in the 
form 

M z 0 y TJ-jE ? Jo /f( 



For a given <p, ft is given implicitly as a function off 

in (29). Transforming the integral in Ij, so that 77 is 

the variable of integration, the quadrature can be 
analytically carried out and we readily obtain 



J2 \ -z 0 '<p W 



sm 



This equation is easily inverted, so that, using the 
preceding expressions, we obtain iJfi = Hy) or u i =f(y) 
explicitly, as has been done in (33) below. 

In collecting the most important expressions for the 
final, complete solution, it is desirable to include the 
expression for the so-called "boundary-layer Reynolds 
Number" B s . If Ui is the actual (dimensional) 
potential velocity just outside of the boundary layer, 




Figure 5.— Alternative picture of the method of joining inner and outer solutions. 

and 5i is the corresponding actual boundary-layer 
thickness, then Us is defined by 



5 5 = 



In presenting the final results it is also convenient to 
use a slightly different picture from that discussed in 
connection with figure 4. For conceptual purposes 
the displaced origin for the outer solution seems most 
suitable to the authors, since it gives a continuous final 
solution in the f (£) or f (^) plane. The notation for 
the corresponding expressions giving results in the 
u(y) plane is, however, a little complicated. The 
alternative picture indicated in figure 5 gives identical 
results to that of figure 4 but leads to simpler expres- 
sions for u and y. In figure 5 we have made the 
origins for f < and f 0 coincide so that we have only one 
£ variable instead of the two used formerly (£ t and 
£ M ). At the joining point tj> there is now a discon- 
tinuity of magnitude e between the values of £ for the 
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inner and outer solutions (£ w — kaj=e). The complete 
solution is obtained by following the inner solution 
curve from the origin to ^ and the outer from £, to z Q . 
(cf. the solid curve.) The solution u(y) is then con- 
tinuous although f(S) is discontinuous, as is clear from 
the formulas given just below. 

The essential results of the entire foregoing analysis 
may be summarized as follows: For any given value 
of a between the stagnation point origin and the 
separation point: 

Writing y * = y^/B^z?, g= -yj tb 
for f i.e., uHuj'. 



then 



^=|| 1-VT+^sin 



V u 

for ^tu i-e., v&*hi ^=^tt) and y* 

y=S, i.e., y* =8* for £=f s where 
^(1= Is) = 0.99 



(33) 



At the separation point, 5=0, we have the following 
simplifications: 



„ u Ujll — COS ~r= | 



(33a) 



It will be noticed that the boundary-layer thickness 
y =5 has been denned in accordance with (21) as that 

value of y for which ^ =0.99, i.e., -^=0.995. This 

numerical value was chosen because of the excellent 
agreement it gave with the Blasius expression for 5 in 
the case of flow with constant pressure. It will also 
be noticed that the final formulation given above 
completely eliminates the necessity of calculating 
% l} or e. 

Before leaving this general analysis of the boundary- 
layer problem, the method of determining the separa- 
tion point should be indicated. It would, of course, 
be possible to determine the value of <p s giving this 
point by trial and error, using the procedure already 
outlined in this section, until a value of <p was found 

for which B 2 = (^) =0. The problem can, how- 



ever, be much more easily solved directly by the follow- 
ing method. Since B 2 =0, we have from (31) for the 
separation point: 



--1 
2 

V2 



B 2 =0 (34) 



From the outside solution, using the condition that 



where f] and f 2 are known functions of <p. 
employ the conditions: 



We now 




which, substituted into the second equation of (34), give 
— ^ 2 -^=fs(*>), where f 3 is also a known function of <p. 

So 

But zo and z^ are also given functions of <p. We there- 
fore have an equation containing only <p whose solution 
obviously gives <p s , i.e., the value of tp at the separation 
point. In other words, we determine the separation 
point position <p s by solving an equation of the form 

2^ 5 ) = 0. 

The procedure is entirely similar to that discussed 
earlier in connection with the use of the unmodified 
outer solution alone. The calculation of the velocity 
profile, etc., is carried out in exactly the same manner 
as was described above in this section for an arbitrary 
value of <p. 

13. APPLICATION TO THE BOUNDARY LAYER FOR A 
" SINGLE-ROOF PROFILE" 

In this section we consider the application of the 
above analysis to the simplest case of external flow with 
pressure increase given by U 2 =p 0 + Pi<p- We call this 
TPfo) profile a "single-roof profile" for obvious reasons. 
For any actual flow around a solid body we have an 
upstream stagnation point, i.e., ZP^O at <p = 0. In 
order to give our single-roof profile any physical sig- 
nificance we must, therefore, consider it as a limiting 
case of a "double-roof profile" defined by 

E7 2 =6 l¥ > for 0<<pTS.<pi 
U 2 =p a + p x <piov <p>z<pi. 

The limiting single-roof profile is obtained by letting 
«>r^0, while & l¥ >i=W. (cf. fig. 6.) The final 

velocity profile may then be written: 
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and may be considered as a first and crude approxima- 
tion to the flow over the upper surface of an airfoil. 
To fix ideas we consider this velocity profile as such 
an approximation. At the trailing edge of an airfoil 
the velocity has very nearly the value of the undis- 
turbed velocity far from the airfoil. If this be taken 
as the basic velocity used in defining the Eeynolds 
Number of the problem, then we may take U (trailing 
edge) = l in accordance with (3). Once again to fix 
ideas we may consider our airfoil only in the normal 
range of angles of attack, so that the position of the 
forward stagnation point varies very little and may 
be assumed to remain fixed. If we then take the 
distance along the upper surface from stagnation point 



Double roof prof He 




0 <Pi <P 

Figure 0— The single-roof potential velocity profile as the limiting case of a double- 
roof profile. 

to trailing edge as the characteristic length of the 
problem, the dimensionless coordinate x along the 
surface will have values between zero at the stagnation 
point and unity at the trailing edge. From the rela- 
tion between x and <p given in (9), it is then very simple 
to deduce the function U(x) from TP(<p). The two 
alternative expressions giving U are finally 



2z 0 



C7 2 =Z7 l 2 -2(Z7 1 -l)d 

Z7 = E7i -(Ui-l)as. J (6d) 



The corresponding profiles are plotted in figure 7. 

In obtaining the outer solution it appears that the 
results are identical whether we substitute (35) directly 
into (14) or use the double-roof profile with (17) and 
then pass to the limit <pi~^0. In either case, we obtain: 



z»= Ufa® - 2(J7i- l)bb(© +gdQ]v 
= tfi 2 <7o(£)-2(Z7 1 -l)A 1 ®*> 



(36) 



For the separation criterion, on the other hand, we 
must use the double-roof profile and the limiting proc- 
ess. If we carry this out using (18), we easily get 



^""W-i) (37) 

where the subscript w is used to indicate the fact that 
the value of <p s is obtained from the separation criterion 
U 




Trailing 
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Trailing 
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Figure 7.— The final single-roof profile as function of z and <p. 

for the outer solution only. The functions g 0 (£), <7i(£), 
and A-i(£) are plotted in figure 8. 
To carry out the joining procedure with the inner 

solution, we need — = 1 — -> and in the expression for 
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Figure 8— The universal functions go; n, and hi. 



2.5 



this quantity it appears .convenient to replace <p or U 2 
by a new variable denoting position along the surface. 
We, therefore, write 

<7x 2 1 
«= w = 



n U-l - 



' (38) 
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and obtain at once from (35) and (36) 

^=1-2A 1 (9+2« ?1 (|) 
So 



(39) 



d 2 



From the condition -^J^ (if=f Mf ) = 0,wegeta=a(£ u/ ). 

Since is always small, the easiest method which has 
been discovered for carrying out this calculation in 
general is to expand the functions g^) as power series 
about 1=0, using the expansion for the probability 
integral given in section 8. In the present case it was 
found possible to neglect terms of degree higher than 

5 in | in the expression for — > and this is probably also 



true in general. Having a = a (&,,-), it is easy to cal- 




culate 0r). and J J fr° m (39) as functions of 

Turning now to the inner solution, so as to determine 
the separation point, we obtain from (35) and (38): 

(40) 



-Zo> 



Zo 



■a-1 



so that (34) becomes 



Now replacing (^J* by (~). as obtained in the last 



paragraph, we have 



= f(? u ,), since a=a(&, ; ). 



(say), which may be 



calculated and plotted. In accordance with the con- 
siderations presented at the end of the preceding sec- 
tion, the value of for which i^i (&>,-) vanishes, is that 
corresponding to the separation point. The separa- 

tion point values of a and — are given by this value of 

Zo 

£<•>>■• (?'i)s is then determined using (34) and e, from (32). 
The numerical values obtained in this way for the 
"single roof" profile are: 

« s = 1.241, (&) =0.246, (&,,) =0.262, e s = 0.026 (41) 

To get a picture of the complete boundary-layer be- 
havior, we must consider other values of a between 1 
and 1.241, corresponding to points along the surface 
between the stagnation and separation points. For 
such intermediate points the procedure is the following: 
We assume a definite value of a and this, in view of 
the calculations outlined above, immediately gives 



1 expresi 



and I f ' Substituting the last two into 
the second expression of (31), remembering that 



we obtain B 2 . Although they are 

not needed for the determination of boundary-layer 
velocity profiles, we could then calculate £ t) and hence 
e from the first equation of (31). 

The final velocity profiles for any of the values of a 
considered are now determined from (33) or (33a). 
The integral occurring in these expressions will usually 
have to be evaluated by an integraph or some other 
graphical or tabular means. The complete calcula- 
tions have been carried out for a series of three values 
of ex. In order to visualize the results and to present 
them in a form independent of the magnitude of U u 

they are given in terms oijj> V*, and x* where 



Y T n=> y* = y V# V^i-i", = * 



The velocity profiles and boundary-layer thickness in 

terms of these variables are indicated in figure 9 by the 

solid curves. The joining points are shown by the oval 

symbols with the j. It is very interesting that in the 

u, y diagram the inner solution is so large a fraction of 

the complete solution, while in a f , £ plane the outer 

solution goes nearly to the wall. It is just the small 

values of $ at the joining point which make it seem that 

the approximation of our joined solution to an exact 

solution of the boundary-layer equations is satisfactory. 

For the sake of comparison the analogous calculations 

have also been carried out using the outer solution only. 

The results are given by the dashed curves in figure 9. 

At the separation point as determined from the outer 

solution (cf. equation (37) the boundary-layer thick- 

>- 

ness becomes logarithmically infinite since — is pro- 

Zo 

portional to £ 2 near £=0. This seems to be a general 
result and indicates the unsuitability of the outer 
solution alone for the discussion of separation-point 
characteristics. On the other hand, the outer solution 
is very close to the joined solution for the first profile 
(going downstream from the stagnation point), and is 
not very different for the second profile. It would 
appear, therefore, that the outer solution alone prob- 
ably represents a fairly satisfactory approximation in 
general, as long as the immediate neighborhood of a 
separation point is not approached. 

For the single-roof velocity distribution discussed in 
this section, the Pohlhausen method of approximate 
integration of the boundary-layer equations also leads 
to simple expressions. In deriving these expressions 
it is again necessary to consider the single-roof profile 
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as a limiting case of a double-roof profile, but this 
limiting process is very easily carried out, and leads to 
the following simple form of the Pohlhausen equations: 



the outer solution is known to be inexact near a sep- 
aration point, this agreement serves to cast some doubt 
on the Pohlhausen solution in such a region. 




.04 . .06 .08 J0< .13 .14 '.IB 29') 30 

Separation point.-'' Separation point. Separation pointy 

complete solution Pohlhausen solution outer solution 

Figure 9.— Velocity profiles and boundary-layer thickness for the single-root potential velocity distribution. 



logo <7 1 -(S-l) a; =loge ^ =Ix 



_ (*»_. 



-213.12 + 5.767? + ^ 



■ 7257.6 + 1336.3277- 37.92tj 2 - 0.8ij 3 dl} 

i2+\ /y*\ _ wv _ 4-x/j/*v e- yy*Y 

6 \S*J 2\S*J 2 \8*J + 6 \S*J 
5* = 



(43) 




Figure 10— The double-roof profile. 

The function Ix(X) may be evaluated analytically or 
more simply with an integraph. For the separation 
point X = — 12 which gives a s = 1.404. Boundary-layer 
profiles have been calculated from (43) for this value 
of a as well as for the two values considered above 
with our joined and outer solutions. The results are 
plotted in figure 9 as dotted curves. The agreement 
with our solutions is quite good for the first two pro- 
files, and at the third profile the outer and Pohlhausen 
solutions are in remarkably good agreement. Since 



The expressions for Bs s as calculated by the three 
methods are: 



Rs s (joined solution) = 2.18 • 



'Vtfi— 1 

Bs 5 (outer solution) = co 

B h (Pohlhausen) =2.93 -7=== V# 



(44) 



14. THE BOUNDARY LAYER FOR A "DOUBLE-ROOF 
PROFILE " 

The results of the last section are particularly sim- 
ple but the single-roof profile is too restricted to serve 
as a satisfactory approximation to many potential 
velocity distributions. The double-roof profile, in- 
volving one more free parameter, is considerably more 
flexible, and may have a much wider field of applica- 
bility. In this section, therefore, the separation-point 
characteristics of a double-roof profile are investigated. 

The notation adopted is defined in the following 
equations, and illustrated graphically in figure 10. 



For 0<<p<<p u U 2 =b<p 

For IP = (6 + 0) n - ftp 



(45) 



We are primarily interested in separation phenom- 
ena and will therefore restrict our considerations to 
the region ¥>>pi, in which region only can separation 
occur. From equations (16) and (19) we obtain for 
the outer solution corresponding to the flow of figure 
10: 

*„ (ft *) = { (& + 0) <d ~ Mh* + b 9 Qh - hi*) , 
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where hi=h/-j=)i hi*=h/ <— — \ and Ai(a:) is de- 
fined in (14), (19), and plotted in figure 8. It is 
convenient to replace z a by — in the following, where 

Zo 

We get at once &j = 1 - 2^*+^ At*)} (46) 



As above, Ave write -7==? 

vV 



so 



that ^=^(1), A 1 * = A 1 (i^_2_) 



Considering — as a function of ? and restricting our 
attention for the moment to variation with p, it ap- 
pears that there are in (46) two variables: and 

— - — • However, the relation which must exist be- 
<p—<Pi 

tween these two variables is very easy to deduce. For 

TP (6 + 18) ft- fr» <Pi P v-Vi 
b<p b<p <p b <p 

<p 0 



or 



b<p 

Equation (46) may therefore be rewritten: 
f u = l-2{A 1 *+^(A 1 -A 1 *)} 



hi = h^), hi 



where #=^p> r =|> and % = ~J= 




(47) 



The variable &, defining position along the <p axis, 
and the parameter r, giving the ratio of the slopes of 
the two portions of the double-foof profile, recur con- 
tinually throughout the succeeding analysis. 

The first problem is now to investigate the joining 
of the above outer solution with the inner solution so 
as to determine the location of the separation point. 
Equation (34) gives, for the inner solution at the 
separation joining point 

52 £i 
where ~§jjr = 0 



smce 



t 0.404/rA* 

^=2.8287^)^ 



(48) 



Lo I _1 go I 



To join the two solutions, we must now require that 

0 and put f u , = f f , = 

In view of the fact that ^ is small, the most satisfac- 
tory method which has been found for proceeding with 
this calculation is to use the power series expanion of 
A x (a;) about x=0. Using the expansion which is given 
in section 8 for the probability integral occurring in 
h u the espansion is found to be 

h(x) = 0.5000- 1.1284 x+a; 2 -0.3761 a?+ 

The next term in the series is of degree 5 in a; and its 
neglect has been found to cause no perceptible errors 
in the results, so that for purposes of calculation the 
series is cut off after the term in x*. Substituting this 



expression for hi(x) in (47) and putting 



m S I _JLo I = 0, 



wo 



obtain % aj =fi(&, 7*). This in turn leads to expressions 
of the form 

(£2) =f 2 (#, r) and ( ) =f 3 (#, r) 

At the separation point tf = # s . Using these relations 
and the second equation of (48), the conditions 



f«i=f <; ,=fr and 



lead to a relation of the form & s = # s (r). 

From equations (47) and (48) we now determine the 

separation point values of and e = £ <y — £ UJ as 

functions of the parameter r. The numerical results 
are indicated by the curves of figure 11. 

Before proceeding to the determination of the actual 
position of the separation point, let us consider for a 
moment the value of Rs s — the boundary layer 
Reynolds Number at the separation point. From 
equations (33) and (33a) this is given by 



+ J2r&. 



^(f=| 5 )=0.99. 

So 



(49) 
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In carrying out the numerical calculations the integral 
in 5„* was evaluated with a Coradi integraph for a 
series of values of # s , i.e., a series of values of r. 

The essential results of the analysis to this point may 
now be presented in the form # s = #,(r) and 5*=S*(r); 
the results are therefore expressed in terms of the 
single parameter r. In proceeding to the final results 
(separation-point location and value of Bs s ) it is 
necessary to introduce one more parameter which we 
take to be <pi. From the original relation U 2 = (6 + 0) <p x 
— p<p, and from the definition of & it follows immedi- 
ately that 

1+r „ Z7 2 1 + r 
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I'ioure 11. — Soparation, joining-point characteristics lor tho double-roof profile 

If we therefore write 



1 + r 



Xs Vr(l 



2(1 +r) 



5* 



(50) 



where 0 S and A 5 * are functions of r only, then 

The functions B s and A s * are plotted in figure 12 as 
functions of r=|- Figure 12 and equations (51) give the 
position of the separation point (<p s ) and the corre- 
sponding value of -j= for any double-roof profile, in 



a 

terms of the two parameters r= ^-and <pi. The extreme 

simphcity with which all of the major separation-point 
characteristics may be determined, for the rather 
general class of potential flows considered, is quite re- 
markable. 1 

Two rather interesting results of the above analysis 
are worthy of special mention. First the small magni- 
tude of the dependence of 8 S * on r is illustrated by the 
curve of 8 S * vs. r which is included in figure 12. In 
view of the definition of 5* from equation (33) 

where S s is the actual dimensionless boundary-layer 
thickness, this nearly constant value of 5 S *=2.5 
means that we have 
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Figure 12.— G oneral lunctions giving soparation-point characteristics lor the double- 
roof profile. 

In other words, for a given Reynolds Number B, the 
thickness of the boundary layer at the separation point 
is approximately inversely proportional to the slope of 
the second portion of the TP(<p) (or U(x)) profile, and 
independent of the maximum value of 77 or of the man- 
ner in which it is reached, i.e., independent of 6 or tp t . 
The second interesting result lies in the type of simi- 

■ Since 0. behaves approximately like y as r approaches zero, it would appear that 
the function rS. might have some physical interest. Actually it is easy to show that 
^£=1— r (<?,— 1) and this function does have an obvious physical significance. For 
the limiting case r=0, which corresponds to the single-roof profile, we have found 



I^=i=0.806 (cf. equation 41). 

Ul* ott 



As r increases increases practically linearly to 

a value of 0.SS0 at r=0.25. In other words, as the distance along the boundary from 
the stagnation point to the point of maximum potential velocity increases, the po- 
tential velocity at which separation occurs, approaches more and more closely to the 
velocity. This is exactly the behavior which simple physical considera 
Hons would suggest. 
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larity law which is found to hold for a family of double- 
roof profiles. For, if we consider a family of such 
profiles, all of which have the same values of <pi and of 

r as in figure 13, then the values of <p s and of —j= are 

identical for all the members of the family and for all 
values of R. This may be considered as an extension 
of the well-known general boundary-layer theorem 
which states that for any given potential flow function 
the position of the separation point along the surface x s> 
Rs 

and the value of -j= are both constant and independent 
of R. 




Figure 13.— Family of similar potential velocity distributions giving the same 
values of p. and ot B,J-[R. 

CONCLUSION 

In the above sections we have developed a new 
method for determining approximately the characteris- 
tics of a laminar boundary layer. General formulas 
and procedures have been explicitly given enabling the 
direct application of the method to a rather general 
class of external flow functions. Two examples are 
considered in detail and represent approximations to 
the type of flow which takes place over the upper sur- 
face of an airfoil. In most practical cases the boundary 
layer is turbulent over the main portion of the airfoil 
section, but the characteristics of the small laminar 



region downstream from' the stagnation point have a 
great influence on the "transition phenomena" which 
determine the maximum lift at low and moderate 
values of [Reynolds' Number. In the first example, 
the calculation of the boundary layer starts from the 
point where the potential velocity has its maximum 
value, and the development of the boundary layer and 
the location of the separation point are investigated in 
their dependence on the rate at which the velocity de- 
creases with the distance from the point of maximum 
velocity. In the second example, the calculation starts 
from the stagnation point and the effects of both the 
region of increasing velocity and that of decreasing ve- 
locity are taken into account. In the latter example, 
the numerical calculations are restricted to the charac- 
teristics of the boundary layer at the separation point. 

The authors wish to express their indebtedness to 
Dr. N. B. Moore for his kindness in checking the der- 
ivations of the formulas of this paper. 



Guggenheim Aeronautics Laboratory, 
California Institute op Technology, 
Pasadena, Calif., May 22, 1984. 

REFERENCES 

1. Pohlhausen, K.: Zur Naherungsweisen Integration dor 

Differentialgleichung der laminaren Grenzsohicht. Z.f. 
A.M.M., 1921, p. 252. 

2. Mises, R. v.: Bemerkungen zur Hydrodynamik. Z.f.A.M. 

M., 1927, p. 425. 

3. Hiemenz, K.: Dissertation Gottingen (1911). 
Boltze, E.: Dissertation Gottingen (1908). 

Thorn, A.: The Boundary Layer of the Front Portion of a 

Cylinder. R. & M. No. 1176, British A.R.C., 1928. 
Green, J. J.: The Viscous Layer Associated with a Circular 

Cylinder. R. & M. No. 1313, British A.R.C., 1930. 
Falkner, V. M., and Skan, Sylvia W.: Some Approximate 

Solutions of the Boundary Layer Equations. R. & M. 

No. 1314, British A.R.C., 1930. 
Goldstein, S.: Concerning Some Solutions of the Boundary 

Layer Equations in Hydrodynamics. Proceedings, Cam- 
' bridge Phil. Soc, vol. 26, part I, 1930. 

4. Millikan, C. B., and Klein, A. L.: The Effect of Turbulence. 

An Investigation of Maximum Lift Coefficient and Tur- 
bulence in Wind Tunnels and in Flight. Aircraft Engi- 
neering, August 1933. 

5. Prandtl, L.: Gottinger Nachrichten, 1914, p. 177. 



